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I. INTRODUCTION 

In recent years, experimental advances in the domain of cavity quantum electrodynamics 
(QED) and trapped ions has motivated increasing interest in the engineering of nonclassi- 
cal states. The preparation of quantum states constitute a crucial step towards testing 
fundamentals of quantum mechanics, such as nonlocality and decoherence 0, and to the 
development of devices such as logic gates for implementing quantum computation [0 . The- 
oretical schemes have been put forward for engineering arbitrary states of the radiation field, 
both trapped in high-Q cavities 0-^] and as a travelling wave ^j^, besides electronic and 
vibrational states of trapped ions p|^T3|. 

In references the desired superposition of photon-number state in a cavity field is 

engineered from the vacuum state, photon by photon. In Vogel et a/.'s scheme [^, resonant 
atom-field interactions are required to create a general single-mode resonator, while Parkins 
et al. ^ suggest a method using adiabatic transfer of Zeeman coherence, and Law and 
Eberly |^ have provided an approach without the need to prepare multilevel electronic 
superpositions. In reference the authors, despite starting from a coherent state of the 
cavity field, do not provide a scheme for generating an arbitrary superposition state. In 
the travelling wave domain, Pegg et al. [0] proposed a method to generate an arbitrary 
superposition of a zero- and one-photon field state, based on state truncation of travelling 
optical fields. A similar scheme for generating arbitrary quantum states of a travelling wave 
was proposed by Dakna et al. 

Concerning the engineering of trapped ionic states, a technique for preparing even and 
odd coherent motional states of a trapped ion, via laser excitation of two vibrational tran- 
sition, has been described by de Matos Filho and Vogel |^ . Starting from the vacuum state 
as in Ref. 0, Kneer and Law ||lOl proposed a scheme for generating a general quantum 



entanglement of the electronic and vibrational states of a trapped ion, besides discussing 
the engineering of two vibrational degrees of freedom of a single trapped ion. Drobny et al. 



Tl| also provide a technique for deterministic preparation of two-mode vibrational states. 
A method for creating deterministic electronic Bell states of two trapped ions was reported 
by Solano et al. [|I^], and Moya-Cessa et al. |]T3[ presented a procedure to generate arbitrary 
discrete superpositions of vibrational coherent states. 

The preparation of a variety of nonclassical states has been achieved in cavity QED 



||l| Jl^JT6[| and trapped-ion phenomena ||T7| , p!8| . "Schrodinger cat"-like states of both the 



radiation field trapped in superconducting cavities and quantized motional states of 
trapped ions have been achieved experimentally. The preparation of pure photon number 



states of the radiation field has been recently reported 



The feasibility of engineering trapped ionic states relies on the fact that decoherence of 
quantum superpositions can be made negligible by suppressing spontaneous emission using 
metastable transitions. Coherence of the atomic population survives for many Rabi cycles 
of the Jaynes-Cummings (JC) interaction [^], during which, at sufficiently low pressure, 
collisions with background atoms can also be avoided [|T^]. Here we mention a recent, 
remarkable experimental achievement, in which the decoherence of superposed motional 
states of a single trapped ion was controlled through engineered quantum reservoirs . 



To accomplish this, laser cooling techniques have been considered to generate an effectively 



zero-temperature reservoir as suggested in Ref. |2T[]. Differently from processes in cavity 
QED, where the cavity- damping mechanism is the source of errors and decoherence, in the 
domain of trapped ions it is assumed that the errors arise from fluctuating electrical fields 
of the trap and intensity and phase fiuctuations in the exciting laser pulses [^]. Stochastic 
models have been proposed to deal with this sources of errors p3|^6|, and they are shown to 
be in good qualitative agreement with recent experiments. Recently, instead of a stochastic 
mechanism, Di Fidio and Vogel proposed a model where the observed damping in Rabi 



oscillations is caused by quantum jumps to an auxiliary electronic level. Other less 



important error sources in trapped ions, such as collisions with the background gas, are also 



present p3|,p8|. 



An advantage of engineering trapped ionic states over cavity-field states is that the 
latter is a more demanding process. In fact, a bunch of two- level atoms are required to 
generate a cavity field, each atom being appropriately prepared by a Ramsey zone (in a 
particular superposition state), with its velocity selected before interacting with the cavity 
field, and detected before the passage of the subsequent atom ||29|. On the other hand, the 



preparation of an ionic state is designed by just switching laser pulses on and off, alternately, 
to manipulate the electronic and vibrational states of the trapped ion. However, we show 
here that the errors introduced by the intensity fluctuation in the exciting laser pulses has 
a severe effect on the engineered state, even more dramatic than the errors introduced by 



the cavity- damping mechanism when preparing a cavity-fleld state ^1 . 

In this paper we present a scheme for engineering an arbitrary motional state of a trapped 
ion by the projection synthesis method, which was originally proposed in the travelling-wave 
domain for the measurement of particular properties of the radiation fleld, such as its 
phase or its Q-function This scheme consists in sculpting an arbitrary motional ionic 



state from a coherent motional state previously prepared in the ionic trap. The technique of 
sculpting an arbitrary state from a coherent superposition through the projection synthesis 
technique was previously developed in the cavity QED domain |2^. However, here we take 



advantage of the facility to manipulate trapped ions to make the sculpture process even 
more attractive in terms of experimental implementation. In the cavity QED domain, the 
sculpture technique considerably improves previous schemes for generating an arbitrary 



state of the radiation fleld [0]. Whereas in Ref. atoms are required to generate a fleld 
with a maximum number of photons equal to A^, the sculpture technique requires about 
half of this number, exactly M = int [{N + 1) /2]. This is due to the fact that we begin 
our process from a coherent state previously injected into the cavity, instead of from the 
vacuum state as in p. So, instead of building up the desired state photon by photon, we 
proceed to sculpt an existing coherent state with an appropriate average excitation previously 
calculated. The situation is analogous to the present proposal for creating motional ionic 



states. Instead of requiring steps to generate an arbitrary motional state with a maximum 
number of phonons equal to [0, our technique utilizes just M steps, as in the cavity 
QED context. 

We stress that in the present work we elaborate the engineering process of the vibrational 



ionic state in the realistic presence of noise. Following the reasoning in Ref. we consider 
the noise arising from the intensity fluctuations of the laser pulses used to manipulate the 
electronic and vibrational states of the trapped ion. Thus, after discussing the fundamental 
interactions between the trapped ion and a classical field in Sec. II, we describe the sculpture 
technique for preparing the ionic vibrational state in the ideal case (absence of noise) in Sec. 
III. In Sec. IV, we calculate the fidelity of a sculpted state under the effects of noise. Here, 
a phenomenological-operator approach pO| , pl| will be applied to account for the presence of 
noise in the engineering process. This approach has been developed precisely to compute the 
effects of noise in complex processes such as quantum state engineering and teleportation of 
quantum states; processes where the protocol requires several steps of quantum interactions. 
Through the phenomenological-operator approach, first developed in the domain of cavity 
QED pD| , pT| , the effect of noise is introduced directly in the evolution of the state vector of 
the whole system, instead of turning to the evolution of the density operator as in the ah 
initio methods. In Sec. V we present a technique for optimizing the fidelity of the sculpted 
state which is based on an appropriate choice of the ion-laser interaction parameters in the 
expression of the sculpted state computed in the presence of noise. Also in Sec. V, we 
illustrate the scheme by sculpting a phase state and computing the fidelity-probability rate, 
a cost-benefit estimate for sculpting a desired state, as defined in Sec. III. Finally, in Sec. 
VI we give a summary and draw conclusions. 

II. MODEL 

We consider one single trapped ion of mass m in a one-dimensional harmonic trap whose 
frequency is u. The ion has forbidden transitions between two internal electronic states 
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(excited | t) and ground | |) states, taken as hyperfine sublevels of the ground state), 
separated by frequency Uq and indirectly coupled by interactions with two laser beams, of 
frequencies uJi and 002, in a stimulated Raman- type configuration. As indicated in Fig. 1, the 
laser beams are detuned by A from a third more excited level |r) which, in the stimulated 
Raman-type configuration, is adiabatically eliminated when A is much larger than all of the 
following: the linewidth of level |r), the coupling associated with the | f) ^ |r) and | |) ^ 
|r) transitions, and the detuning 6 = ujq — ujl iy^L = t^i — ^^2) [HjH]- The transition between 
I I) and a fourth level \d), achieved by another laser strongly coupled to the electronic 
ground state, is considered in order to measure the ionic vibrational state by collecting the 
resonance fiuorescence signal, which is the probability of the ion being found in the internal 



state li) p 



The Hamiltonian that describes the effective interaction of the quantized motion of the 
ionic center-of-mass (CM) with its electronic degree of freedom is [0,0], in a frame rotating 
at the "effective laser frequency" ul { h = 1): 

where ip is the phase difference between the two lasers, a+ = | t)(i 1; o"- = | |)(t | and 
are the usual Pauli pseudo-spin operators, a'' (a) is the creation (annihilation) operator of 
vibrational quanta, Q is the effective Rabi frequency of the transition If) ^ ||) and, rj is 



the Lamb-Dicke parameter defined as |18,|35|: 



V = (2) 

Here A A; = {k i — k 2) ■ i ,\k 1(2) | = uji(2)/c, k i { k 2) being the wave vector for laser 1(2), 
and i is the unit vector in the direction of the trap axis . 

Written H in the interaction picture and then expanding the resulting Hamiltonian in 
terms of the Lamb-Dicke parameter we get 



^ mW. 

.m,l=0 



(3) 



Assuming the Lamb-Dicke limit, for which 77 ^ 1, where the ionic CM motion is strongly 
localized with respect to the laser wavelengths, we obtain the simplified Hamiltonian 

where resonance is achieved by tuning the laser frequencies to obtain 6 = —iu {i = m — I). 



Considering the realistic value for the trap frequency z//27r ^ 11.2 MHz |jT5[|, the optical 
rotating wave approximation leads to the Carrier Hamiltonian {i = 0) when tuning the 
effective laser frequency to obtain 5 = 0: 

H, = Q{a+e-'^ + a_e''^). (4) 

This Hamiltonian induces the transition \n, |) < — > \n, f) (where \n) indicates a motional 
Fock state), and is responsible for rotating only the internal electronic levels of the ion wave 
function in accordance with 

e-'^^^ |n, T) = cos (fir) \n, t) - ie'"^ sin (fir) \n, i) , (5a) 
e"^^^" |n, i) = cos (fir) |n, i) - ie"^'^ sin (fir) |n, ]) . (5b) 

When tuning the effective laser frequency to obtain 6 = —u, the optical rotating wave 
approximation leads to the Jaynes-Cummings like Hamiltonian {£ = 1) corresponding to the 
first red sideband, 

Hjc = in7]{aa+e-''^ - a^a_e'^), (6) 

which induces the transition \n, I) < — > \n — 1,1); iii such a way that the electronic and 
vibrational modes evolve as 

e-'^^^^ \n, T) = a \n, t) - e-'^S^ \n + 1, j) , (7a) 
e-'^'^^ \n, i) = a-i \n, i) + e'^^S^-i \n - 1, T) , (7b) 



where C„ = cos{gTy/n + 1), Sn = sm{gT\/n + 1), r is the duration of the laser pulses, 
and g = fir/. Finally, we note that it is possible to obtain the Anti- Jaynes-Cummings 



Hamiltonian {C = —1) corresponding to the first blue sideband. However, for the purposes 
of the present paper we do not use this specific interaction which induces the transition 
i) < > \n + 1, t). 

III. SCULPTURE SCHEME (IDEAL CASE) 

In this section we show how to transpose the projection synthesis technique from its 



original travelling wave domain to the context of the ionic trap (similarly to Ref. pS 
where the projection synthesis was applied to cavity QED phenomena). We assume the 
trapped ion to be initially prepared with its CM in a coherent motional state and in the 
electronic excited state |t). Such a state, which can be prepared with techniques available 
nowadays |lT7| , p!8| , reads 



|M;(o)) = |a)®|T) = 5^AW|n)®|T), (8) 

n=0 

with Ai°^ = exp(— |a| /2)a"'/ \fn}.. As the desired state is generated from a previous co- 
herent state, we have denominated the present scheme a quantum state sculpture process 
(as originally done in ||29|). In fact, as shown below, our strategy consists of modelling the 
Wigner function of the desired state, through appropriate laser pulses, from that of the pre- 
viously prepared motional coherent state in (H). The carrier (C) and the Jaynes-Cummings 
(JC) laser pulses work as "quantum chisels" on the initial coherent distribution, as shown 
in the quantum algorithm notation depicted in Fig. 2. In the domain of cavity QED, 
the quantum-chisels are played by two-level Rydberg atoms which are made to interact 
resonantly with a coherent state initially prepared in a high-Q cavity p9[ . 

The whole operation is accomplished in three-step cycles, requiring successively: (i) a 
carrier pulse Ci, to prepare the electronic state in a suitable superposition, (ii) a first red 
sideband pulse JC, to entangle the ionic motional and electronic states, and, finally, (iii) 
a sequence of a carrier pulse C2 and a fluorescence measurement of the ionic electronic 
state. The third step constitutes the projection synthesis, which enables us to synthesize 
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the measurement of a particular superposition of the electronic state and, consequently, to 
synthesize the projection of the motional degree of freedom to the desired sculpted state. 
Thus, the projection synthesis technique is here applied to both of the entangled degrees 
of freedom of the trapped ion, electronic and motional (analogously, in the cavity QED 
domain we synthesized the simultaneous projection of the atomic and the cavity field states 
p9|). A fluorescent signal projects the ion into state ||), while the absence of fluorescence 
projects it into state |t). For the present purpose the absence of fluorescence is crucial in 
preventing the occurrence of recoil of the ionic CM motion. So, at the end of each cycle the 
detection of the absence of fluorescence is required for the successful accomplishment of the 
engineering process. We emphasize that the duration of the fluorescence measurement is an 
order of magnitude smaller than that of the JC pulse and about half that of a C pulse in 
experiments involving the usual parameter values [0. For a JC pulse, gr = 7r/2 and for 
a C pulse, f2r = 7r/2, the respective durations being around 2/is and 0.5/iS, while the time 
for a fluorescence signal is around 0.2yus [0. We repeat this three-step cycle M times to 
synthesize an arbitrary desired state |\E'd) = X]^=o'^" 1^) ® It) (^^^ Fig. 2), where it should 
be noted that the electronic state |t) factorizes. The parameters M and A^^^^ (the maximal 
excitation number of the desired state) are related as shown below. 

Considering the kth. cycle of the sculpture process, let us start with the assumption that 
the ionic state after the {k — l)th cycle is 

oo 

|v[;(^-l)) = ^A(/=-l)|^,|). (9) 

n=0 

As the first step of the fcth cycle we prepare, with the carrier pulse Ci, the electronic state 
in the superposition Mp^ (|t) + I3k where Mp^ = (l + ^ and (3k is a complex 

parameter adjusted by pulse Ci, fixed at an appropriate duration and phase of the carrier 
pulse [as in Eqs. (|^,^) ]. In the second step, the first red sideband pulse JC entangles the 
ionic motional and electronic states as follows: 

oo 

|^^'^> = A/'a E^- "'^ (^-'^ 1^' t) - \n + 1, i) + 

n=0 
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|n, i) + e^^'^/J.^fJ, \n - 1, T)) , (10) 

where Cm^ = cos{gTk\/ m + 1), = sin{gTk^/m + 1), and are the A;th JC pulse 
duration and phase, respectively. Next, in the third step, we have to synthesize the projection 
of state (|10D into a particular electronic superposition state (of the kth cycle) = 
■^Efc (IT) + li))> with A^^ = (l + l^fel ) . The complex number Sk results from the fcth 
rotation of the electronic states induced by the second carrier pulse C2. As a consequence of 
the absence of fluorescent signal (which is the case of a successful measurement), the ionic 
state after the projection synthesis is given by 

00 

|^(^^>=A4|T)(x(^)|V^(^)> = X:Ai^Mn,T)- (11) 

n=0 

The coefficients A^'^'' result from those in equation (^) using the recurrence formula aI^"^ = 
N'kTn \ where 



ri^) = [Cl^^ + eA&:i,) K^'^ + e^-^ASi^^Aiy) - e~^-^e, (1 - 5„,o) ^^A^i^ (12) 
and the normalization constant Mk is given by 



A4 



.n=0 



-1/2 



(13) 



The pulse C2 and the absence of fluorescence signal (detecting the state ||)) are needed 
in order to adjust the measurement of the special superposition jx*''^'')? during which they 
play the role of electronic state "polarizers" and "analyzers" , respectively - by analogy with 
light polarization measurement - and allow us to analyze an arbitrary superposition of |t) 



and II). Such a measurement |Q works as follows: The carrier pulse C2 is appropriately 
adjusted so that the electronic state of the ion in the superposition undergoes a 

unitary transformation to the state | ]). Otherwise, the electronic state evolves to | |) if it 
was initially orthogonal to jx*"*''')- After the carrier pulse C2 a fluorescence measurement is 
required to project |t) or ||). In general, the electronic state after the JC pulse will be a 
superposition of these orthogonal states, so that in unsuccessful cases we flnd a fluorescence 
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signal (detecting ||)), and in successful cases we do not (detecting ||)). Once we find the 
absence of fluorescence the electronic state in the kth cycle has been projected on the required 
superposition Ix^*^"*)- Thus, the projection synthesis technique is able to measure observables 
of the form [x*-*^^) which represent projection operators with measurable eigenvalues 

associated with absence (detecting or presence (detecting ||)) of a fluorescent signal. 

Here we stress an important difference between the present sculpture process and that 
used to prepare an arbitrary state in cavity QED |2^. In the latter, the projection synthesis 



is achieved by measuring the required two-level Rydberg atoms (the quantum chisels for 
sculpting the field state), each in a particular superposition state, through a classical field 
and ionization chamber detectors. Hence, after the projection synthesis, the measured state 
of the kth atom is discarded, since it turns to be useless for the cavity QED process, whereas 
the electronic state of the ion, factored as |t) in a successful event, is ready for the next 
cycle. Therefore, the motional state of the ion is sculpted by means of its own electronic 
states. 

Since the projection of a particular electronic state in the Mth cycle is supposed to finish 
the sculpture process, the following equality must be satisfied 

oo Na 

Mm IT) {x'^'W^) = In, T) = In, T) , (14) 

ra=0 ra=0 

requiring that ^ when n > N^+l. This approximation results in a non-unity fidelity for 
the sculpted state. Usually, the fidelity of a given quantum process is defined to account for 



the inevitable errors introduced by the environment due to dissipative mechanisms ||30| , |3T 



However, as mentioned in the Introduction, in the domain of trapped ions it is assumed 



that the errors arise from noise due to fluctuations in the trap and laser parameters |22 



and to treat this error source we consider the stochastic model proposed in |^, where just 
the intensity fluctuations of the laser pulses are considered. In this section we focus on the 
ideal case in which the fidelity is defined to account only for the errors introduced by the 
approximation (A^^ ^ for n > + 1) inherent in this sculpture scheme. So, this fidelity, 
which does not account for the errors introduced by the environment or fluctuations, reads 
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Z^n=0 "n-'- n 



^^|(vl..|^("))f = " ■ (15) 



1=0 



p{M) 



The total probability of successfully sculpting the desired state is "P = Hfeli ^k, where Pk 
is the probability of synthesizing a particular electronic superposition, from 
the kth entanglement between motional and electronic states, Eq. (p!oD. In other words, Pk 
refers to the probability of measuring absence of fluorescence after the second carrier pulse 
in the kth cycle. This probability is given by 



= ATlATlY^Kf ■ (16) 



n=0 



2 

For an appropriate choice of average excitation of the coherent motional state, |a| =na, 
we see from the recurrence formula ([T^) and the definition of coefficients An\ that the 
coefficients T^n ^^ depend on powers of a, varying as / a/jI, with n—M < j < n+M. In fact, 
it is straightforward to conclude that, after one application of formula (|T2D the coefficients 
ri*^^ are proportional to | aI,'^]^ , aI*^ , A^^^^ | ; after two applications of the formula (|T2|) 
it follows that vi^^^ oc | A,[j^2 ; A^^^ A^*^ ^■^ , A^^^j^ , A^^2 | after M apphcations 
we finally obtain T^^^^ oc |a^°2^^, A^°2j^^_^^, Ai°\ A^°|j^^_p A^l^jj^^j. Therefore, from the 
definition of A^'''* (initial coherent motional state), we note that F^*^^ depends on powers of a, 
as mentioned above, and the choice of the average excitation Ua that ensures P{Nii—M+l) = 
I {Nd — M + l\a)f ^ (satisfying the requirement that Ai^'^^ ^ when n > Nd + 1) results 
in a higher fidelity JF at the expense of a lower probability V. In fact, it is evident from the 
denominator of Eq. (p!5| ) that the lower the number / in the sum of significant coefficients, 
the higher the fidelity. On the other hand, in Eq. (0) we observe that the probability V 
is directly proportional to powers of a. As a consequence, Eqs. (|15|) and (|1^) furnish a 
fidelity-probability rate, 71 = J^^V^, a cost-benefit estimate for sculpting the desired state, 
where the parameters ^ and ( are appropriately chosen to weight the contributions of the 
fidelity and the probability in accordance with the aims of the sculptor. In fact, the sculptor 
may decide to privilege the fidelity or the probability in the cost-benefit estimative, and in 
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the present work we have decided to privilege the fidehty, choosing the values ^ = 4 and 
C = 1/2. In order to maximize the rate TZ, we have to play with all the parameters: the 
durations and phases of the JC laser pulses, and ipk, and the carrier pulse (Ci and C2) 
parameters (3k and Sk- We note that a good strategy to maximize TZ consists in starting 
with a choice of rTo, so that P{Nd — M + 1) ^0, and then proceed to maximize the rate 
TZ, increasing ria at the expense of the fidelity. Next, the duration and phase of the JC 
laser pulses and (pk are chosen so as to maximize the rate TZ. Finally, the choice of the 
carrier pulse parameters, (3k and Sk, follows from a particular solution of the equality (p!4D 
which, togheter with the requirement aI^^'* ^ when n > A^^^ + 1, results in the set of A^^^ + 1 
equations 



d 



M 



^Na^ + ^m(3mC)^'^Li ) A^^ 



(A/-1) 



(A/) «(A/-1) 



d^ = Mm \{c^:'^ + .m/3a.«) A^-) + e-^-/3M5i^^)Ai^i-^) - e-'-.M^HAt^r^)] , (17) 



dn=Ui 



M 



To solve the set of equations (|TTp, we apply the recurrence formula (^) M — 1 times in 
order to express the unknown coefficients An^ in terms of the known values of the coefficients 
of the coherent motional state, An\ In this way we obtain a nonlinear system whose free 
parameters are Pi, (3m and Si, ...,eM, which indicate the particular rotation the electronic 
state must undergo in each cycle, during each carrier pulse (Ci and C2). These variables are 
obtained from the known coefficients dn and Ai^"*. The solvability of a nonlinear system can 
be ensured if the number of equations is equal to the number of variables, the latter being 
the parameters of both carrier pulses. One of the equations in system (|1^) has to be used 
to obtain the normalization constant Ilfcli-^fc ^'^d each cycle carries two free parameters 
{(3k, Sk)- Therefore, the minimum number of cycles necessary to guarantee the solution of 
system (^) must be M =int[(A^c; + 1) /2]. This conclusion follows from the fact that in 
our scheme we start from a coherent motional state. The real variables ti,T2, ...,tm and 
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0i,02---0M are used to improve the cost-benefit rate R. 

With this technique, therefore, it is possible to sculpt an arbitrary vibrational state 
by remodelling another initial vibrational state. In particular, we have started from the 



coherent state since it is easily generated [18 



A. Sculpting a truncated phase state 

To illustrate the sculpture technique we now proceed to engineer the truncated phase 
state {Nd = 2) 

2 



* n=0 



(18) 



As mentioned above, the sculpture process for the state I^E^^) requires just M = 1 cycle, and 
we obtain from equations ([T7|) the system 



(19a) 











A(0) + e^^i/5,5«Af)" 





(19b) 



Solving the above system we obtain a fourth-order polynomial equation in the variable ei 
(jSi) by isolating the variable (3i (si) from one equation and substituting it into the other. In 
this way we obtain the roots of the system ( [19a|J19b|) for any fixed set of parameters rT^, gri 
and ifi. When considering more than one cycle to sculpt a state where A^^ > 2, instead of 
Eqs. ( |19a|J19b|) we obtain a set of coupled equations, permitting only numerical solutions 



Following the strategy mentioned above, we start with the average excitation TTq, leading 
to the highest fidelity resulting from P(3) = |(2|a)|^ ~ 0. For our purposes we begin with 
the average excitation TTq, = 0.04. For each value of Ua (choosing a as a real parameter) 
we proceed to calculate the duration (gri) and phase (v^i) of the JC laser pulse which 
maximize the rate TZ. As discussed above, the maximum value of TZ depends on the choice 
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of the parameters ^ = 4 and C, = 1/2^ weighting JF and P, respectively. In Table I we show 
the rate IZ associated with each value of tIq,, from that which maximizes the fidelity (0.04) 
to values that exhibit a continuous decrease of the rate TZ. In addition, four roots {ei,(3i) 
of Eqs. ( |19a| , |19b| ) result when a given pair of parameters [gri^ipi) are fixed, and we have 
to choose the one which maximizes TZ. The value ria = 0.25 results in the highest rate 
TZ = 0.60, which follows from a fidelity JF = 0.99 and probability V = 0.38 of successfully 
sculpting the desired state. 

We display in Figs. 3(a), 3(b) the sculpture process of the desired truncated phase state 
(|T8|) from the Wigner distribution function of the initial coherent state (associated to ria = 
0.25) given by a Gaussian shifted from the origin as W{p, q) = {2/tt) exp [— {q + a)^ — p^] , 
shown in Fig. 3(a). The state (associated with the best rate 7Z = 0.60) obtained after 
one cycle is displayed in Fig.3(b), using the parameters {gTi,ipi,ei and Pi) associated with 
this rate. As we have stressed above, it is possible by the present scheme to sculpt the 
desired state with a higher probability of sucess but at the expense of a smaller fidelity. 
So, the sculpture technique can be evaluated by a cost-benefit estimate, here defined as the 
fidelity-probability rate, which it is up to the sculptor to choose. 

The next section deals with the effects of noise on the sculpture process, i.e., the influence 
of the errors arising from the intensity fluctuations in the exciting laser pulses. In this paper 
we do not consider either the (weaker) effects of the phase fluctuations of the laser pulses 
||23| or the (close to unity) efficiency of detection of the internal state required in the third 
step of the cycle p2|. When the effects of noise on the sculpture process are taken into 
account, following the phenomenological-operator approach, the fidelity defined by Eq. (plSj) 
remains exactly the same, whereas the sculpt field state, after the last required cycle, will 
be entangled with auxiliary states in which noise-operators are defined. So, to the best of 
our knowledge the sculpted motional state will be represented by a statistical mixture pion, 
whereby the fidelity turns out to be JF = [ip^ \pion \ i^d)- 
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IV. EFFECTS OF NOISE ON THE PROCESS 



As mentioned above, we consider here the noise arising from the intensity fluctuations in 
the exciting laser pulses (carrier and Jaynes-Cummings) required to manipulate the internal 



and external states of the trapped ion p3| , |28[| . We will also propose an alternative way of 
engineering an ionic motional state in the presence of noise, which consists in maximizing 
the fidelity of the experimentally achieved state. A phenomenological-operator approach, 
originally developed in the context of cavity QED pn| , PT| , will be considered here to account 
for the evolution of the ionic states under the influence of the fluctuating-intensity laser 
pulses. The strategy is to provide a straightforward technique to incorporate the main results 
obtained by standard ab-initio methods for treating errors in trapped ions (in particular we 



consider the master equation calculations in Ref. ||2^). In short, we introduce an auxiliary 
state space where noise operators are defined to allow for the effects of noise explicitly in the 
evolution of the state vector of the whole system comprehended by the ionic and auxiliary 
states. After computing the evolved state of the whole system, the reduced density matrix 
of the ionic system can immediately be obtained by tracing out the auxiliary variables. The 
phenomenological approach is constructed so that the reduced density matrix turns out to 
be exactly the same as the one obtained by standard methods. It is interesting to note that 
the phenomenological approach resembles the Monte Carlo wave function method in 
the sense that we work directly with the wave function, providing an efficient computational 
tool. 

First we give a short description of the master equation treatment by Schneider and 



Milburn of the laser fluctuations as a stochastic process, the rise or fall of the laser 
intensity being defined as a real Wiener process. After the pulses required to generate the 
desired state, the noise introduced is effectively averaged and the master equation describing 



the ionic system follows from the stochastic Liouville-von Neumann equation ||23 

d 



dt 



p{t) = -i n,pit) 



r 

2 



(20) 
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where Ti. is the interaction Hamiltonian for the C (^) or JC (|^) pulse and the parameter F, 
to be obtained phenomenologically, scales the noise. Considering the subspace composed of 
the eigenstates of Ti, it is straightforward to solve Eq. (PO]), obtaining 

{^i\p{t)\^t) =exppt($^-$±)-2n($^-$;^)' X 

x(<f^|p(0)|$^> (21) 

where 1"^*^) ^^e eigenstates and are the eigenvalues for the interaction Hamiltonian 
described by Eqs. (|D and (|^). For the C Hamiltonian, 

|$^> =-^(|ni)±e^^|nT)), 
K = 

while for the JC Hamiltonian, 

|$^> =i=(|n,i)±ze^<^|n-l,T)), 

Employing the reasoning of the phenomenological-operator approach, we have introduced 
the noise-operators C and J', to account for the intensity fluctuations in the C and the JC 
pulses, respectively. These operators are supposed to act on the auxiliary states |C) and | J), 
so as to the effects of noise explicitly in the evolution of the whole system, now composed 
of the ionic and auxiliary states. After the ion-laser interactions required for the sculpture 
process, the reduced density operator of the ion is obtained by tracing out the auxiliary 
space describing the noise-sources. 

With the noise-operators defined in the auxiliary spaces, we observe that the coupling 
of a general ionic state to the fluctuating C laser pulses evolves in time as follows: 

J2 1^ i) + 1^ t)) ® |C) \^c) = {'^nCn,il {t, ^) + PnCn,n (t, ^)) |n i) ® |C) + 

n n 

[anCn,i^ (t, ip) + l3nCn,n {t, v)) \n ^) ® |C) (22) 

while for the JC pulse it follows 
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J2 1^' i) + 1^' T)) ® I J) ^ I^JC) = J] [«n^n,U (t, 1^, i) + l^nJnM {t, if) \n + 1, j) + 

n n 

anJnU {t, \n - 1, T) + (3nJn,n (t, ^) \n, T)l ® I J) . (23) 



Following the reasoning of the phenomenological approach, the matrix elements 



C 



C ) and ( J 



Jnjk (t, ^) JIj,,, (t, if) I j) (n, n' = 0, 1, 2...., J, f =T 
, I and k, k' =], |), which result from tracing the density operator associated with Eqs. (p2D 
and (p3D over the auxiliary spaces C and J (Trc(jc) \^c{JC)) {^c{JC)\ = Pc{JC))i 
ferred from the standard stochastic method in [|23|. Comparing the reduced density operator 
Pcfjc) with those emerging from the evolution in Eq. (|2T|) of the density matrix associated 
with the state (a„ |^ i) + /5n |^ T))) we obtain the required matrix elements shown in 
Appendix A. Such matrix elements can now be directly applied to any process that involves 
interactions of an ion with fluctuating intensity laser pulses, removing the necessity to per- 
form the typically extensive ah initio calculations. In the case of a process requiring several 
laser pulses (as when sculpting a state with a large Nd) it turns to be practically prohibitive 
to compute the evolution of the ionic system by standard techniques. 



A. Estimating the errors introduced by C and JC pulses 

Let us now consider the initial state |'?/'c' (0)) = |'?/'jc (0)) = |n, |), which evolves, in 
the ideal case (without fluctuations), following Eqs. ( ^b]) and (^) for carrier and Jaynes- 
Cummings pulses, respectively: 



\^c{t)) 



-iHc 



- (0)) = COS (fit) \n, i) - ie-'"^ sin (fit) |n, T) 



l^jc (t)) = e-'^'^^ l^jc (0)) = cos {gtV^ \n, j) + e*'^ sin {gtV^) \n - 1, T) , 

where t and ip are the duration and phase for both pulses. Now, given realistic intensity fluc- 
tuations in the laser pulses, the evolution of states {ipc (0)) and \ipjc (0)) is easily computed 
with the phenomenological approach, leading to the results 

i^c (t)) ® |C) = [Cn,ii (t, if) \n, i) + Cn,i^ (t, if) \n, T)] ® |C) , 



i^jc (t)) ® I J) = [jn,il (i, ^) \n, i) + Jn,u ^) 1^ - 1> T)] ® I J) • 

The reduced density operator of the ion (obtained by tracing out the auxihary spaces C and 
J) are 



%it) =Trc|^c(t))|C)(C|(^cW 



Cn,ii {t, if ) Cl^^^ it, y.))^ |n, i) {n,l\ + (C„,u it. V))^ \n, [) (n, T| + 



+ 



p^^(t) =Trj ^jc(t))|J)(J|(V^c(t) 

^nat {t. v) Jill (t, ^))_^ \n - 1, T) il + (i„aT V^) ^i,iT v))^\n-i, T) - 
From the matrix elements shown in Appendix A and after a straightforward calculation 



we obtain the following fidelities of the evolved states 
to the ideal evolution of \ipc (0)) and I'lpjc (0)) 

0'c=Uc{t) %{t) i^cit) 



ipc (t) ) and ijjjc (t) ) with respect 



1 1 
- + -e 

2 2 



JC 



^JC [t) Pjc (t) ^JC (t) 



1 1 



-2nTg'^t 



for carrier and Jaynes-Cummings pulses, respectively. Note that for F = we have Tc = 
Tjc = 1, as expected. Besides, for n = we have Tjc = 1 since the intensity fluctuations 
in the JC pulse do not exchange energy with the ionic system, maintaining the initial state 
\n = 0, 1) unaffected. As defined in Sec. II, g = t]Q, and the Lamb-Dicke parameter used 
in experiments is rj = 0.202 [jT5[; thus, the noise introduced by the JC pulse has less effect 
than that coming from the C pulse for n < 25. Thus, when sculpting the phase state ( p!8D 
or in any process involving a small number of phonons, such as engineering a motional qubit 
state Co |0) + Ci |1) for quantum computation, as performed in the main source of errors 
will undoubtedly be the C pulses. 
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B. Sculpting a truncated phase state in the presence of noise 



Now we proceed to incorporate the intensity fluctuations of the laser pulses used in 



sculpting the truncated phase state ([ISD which is obtained after one cycle of the sculpture 
process. We assume, for simplicity, that there is no noise in the preparation of the initial 
coherent motional state of the ion. After the preparation of such a motional state and the 
first carrier pulse it follows that: 

oo oo 

^A(°) |n,T) ® |Ci) — J]Af (c;,u (ti) \n i) +C;,n (ti) \n T)) ® |Ci) 



n=0 



n=0 

y 



where ti is the duration of the first carrier pulse and Ci stands for its auxiliary space. To 
simplify the notation we have omitted that the operator C depends on the phase of the 
C laser pulse. Note that in this first step the required superposition of the electronic state 
Nj3 (IT) + h ID) is obtained with a non-unity fidelity due to laser fluctuations. In the second 
step, the first red sideband pulse JC entangles the ionic motional and electronic states as 
follows: 



® I J) ^ 5^ [Xii {t2) CnM (ti) \n, i) + Jnn (^2) {h) \n + 1, i) + 

n=0 

+Jni^ (ts) Cnu ih) |ri - 1, T) + J„TT (^2) (ti) \n, T)] ® | J, Ci) 



where ^2 and v?2 are the duration and phase of the JC laser pulse. Finally, after the second 
C pulse we have 

00 

C2) ^^Af [(c„,u(t3)i^ai(^2)C„,n(ti,¥^i)+C„,n(t3)X^ Ki) 

n=0 

+Cn+l,ll (ts) JnU (^2) Cn,11 (tl) \n + 1, i) + (tg) Xu (^2) (ti) \n + 1 

+Cn-l,n ih) Jni\ (^2) C„,Ti (il) 1^ - 1, T) + Cn_l,t| (ts) ^njT (^2) Cn^U (tl) |n - 1 

IC2, J, 



+ C^at (^3) Jn,u (^2) (tl) + C„,n (ta) X,TT (^2) C„,tT (tl) |n, T) 
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with ts and (f^ standing for the duration and phase of the second C laser pulse. In the 

tp'" ) evolves to 



absence of fluorescence signal, the above entangled state 

oo 

i^") =^Y1 i^nlAl^ ih) Jn^l,n it2) C„„i,tT (tl) (1 - Sn,o) + 
n=0 

ih) + AWC„,it (^3) Jn,U it2) Cn,n ih) + 

+Af C„,tT (^3) Jn,n ih) Cn,n (ti)] 1^ T) ® IC2, J, Ci) , (24) 

and tracing out the auxiliary spaces C2, J and Ci, the ionic reduced density matrix follows: 

00 

Pion(ti,i2,t3) = Tr^2,5,Ci i^'^/ii^'^ = Pn,m{ti,t2,h)\n]) {m]\, (25) 



n,m=0 

where the coefficients pn,m are given in Appendix B and the normalization constant reads 
^ ~ iYl'^=Q Pn,n{ti,t2,t3))~^^'^ . The probability of detecting absence of fluorescence is P = 
and the fidelity of the sculpted mixed state is given by 

:F= {^d\pion{tut2,t3)\^d). (26) 

We stress that in the present situation we do not have a system of nonlinear equations as in 
(p!7|). However, the best parameters (average excitation of initial coherent state, interaction 
times and phases of the pulses) to reach the desired state are determined by numerical 
maximization of the fidelity expression in Eq. (P^D, as shown below. 

V. ALGORITHM FOR OPTIMIZING THE FIDELITY OF A SCULPTED STATE 

We begin by stressing that there is a crucial difference between the sculpture scheme 
for the ideal case, described in Sec. Ill, and for the realistic situation with the effects of 
noise, described in the previous section. In fact, in the protocol proposed in Sec. Ill, it 
is necessary to solve a nonlinear system of equations for jSk and Sk, which arise from the 
recurrence equations relating the amplitudes A^."'^ On the other hand, when taking into 
account the intensity fluctuations in the laser pulses, there is no immediate way to extract 
useful information from a recurrence relation for the operators C„jfc(^i); <Jn,jkit2)^ and 
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Cn,jk (ts) = 0, 1, 2...., j =1, J, and k =t, i) associated with each required cycle. Thus, the 
first question that arises is how can this problem be overcome? We propose the following 
solution: consider a specific desired motional state to be sculpted and calculate the fidelity 
of the sculpted mixed state (^61) . The resulting expression for the fidelity is considerably 
involved and a numerical calculation must be performed to obtain the values tj and ipi 
{i = 1,2,3) that maximize the fidelity-probability rate TZ = J^^.V^. Starting then from 
the average excitation TTq, that leads to the highest fidelity deduced in the ideal case, we 
proceed by numerical optimization to find the value of fia which maximizes the rate IZ. 
Following this procedure, to sculpt the phase state ([T8| ) under the effects of noise, we find 
that fia = 0.25 (coincidently the same value as in the ideal case) results in a highest rate 
TZ = 0.64, which follows from a fidelity J-' = 0.91 and probability V = 0.86 of measuring 
absence of fluorescence. 

Fig. 4(a) displays the Wigner distribution function for the desired state (|T^), which is to 
be compared with the Wigner function in Fig. 3(b) for the sculpted state in the ideal case, 
with fidelity JF = 0.99. Fig. 4(b) displays the mixed state sculpted in the presence of noise 
using the same parameters Pi and ei of Fig.3(b). By comparing both figures we see that the 
interference terms leading to negative contributions in Fig. 4(b) have diffused into each other 
and practically cancel out in Fig. 4(b), which correspond to JF = 0.85 and V = 0.40. Finally, 
Fig. 4(c) displays the mixed state sculpted through our numerical optimization procedure. 
The fidelity and probability, = 0.91 and V = 0.86, respectively, resulting from our 
optimization procedure are considerably higher than those obtained for the mixed state in 
Fig. 4(b). It is worth noting that although the Wigner distribution in Fig.4(b) seems to 
be closely to the distribution in Fig.4(a) than that in Fig.4(c), the optimized mixed state 
possesses a higher fidelity, which is just a measure of the distance between the vectors 
representing the desired and obtained states in Hilbert space and can not be inferred from 
the Wigner function (see discussion in Appendix C). 



To plot Figs. 4 (a)- (c), we have used realistic values [|15| for the Lamb-Dicke parameter, 
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0.202, Rabi frequency, Vt/2Ti ^ 475 kHz, and we assume the parameter F ^ 10 ^s. 



estimated in |23| to obtain a good agreement with the experimental results in Ref |15[]. 



VI. CONCLUSION 

In this work we present a technique to engineer arbitrary ionic motional states employing 
projection synthesis. This method consists in sculpt the desired motional state from a coher- 
ent motional state previously prepared in the ionic trap. The sculpture of arbitrary states 
from a coherent superposition by projection synthesis was previously developed in the cavity 



QED domain p9[. However, here we take advantage of the facilities in manipulating trapped 
ions to make the sculpture process even more attractive for experimental implementation. 
As in the cavity QED context, instead of requiring N laser pulses to generate an arbitrary 
motional state with a maximum number of phonons equal to N [^, our technique utilizes 
just M = int [{N + 1) /2] pulses. 

The sculpture scheme was also developed for the realistic situation of intensity fluctua- 
tions in the required laser pulses. In this connection, a phenomenological-operator approach 
developed recently to account for errors in complex quantum processes was applied [l30| , |3l| . 
This approach furnishes a straightforward technique to estimate the fidelity resulting from 
the engineering process, withoud the need to perform the usual extended ab initio calcula- 
tions required by standard methods. The reasoning behind the phenomenological approach 
developed here is to incorporate in a concise algebraic form the main results obtained from 
the application of the master equation (to the effects of noise arising from the intensity fluc- 
tuations in the laser pulses) as done in Ref. p3[. By defining an auxiliary state space where 
phenomenological operators are defined, we account for the effects of noise explicitly in the 
evolution of the state vector of the whole system comprehended by the ionic and auxiliary 
states. After computing the evolved state of the whole system, we can immediately obtain 
the reduced density matrix of the ionic system by tracing out the auxiliary variables. We 
stress that the results obtained here from the phenomenological procedure are quite general, 
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and in principle can be applied to any quantum process in an ionic trap, such as quantum 
communication, logic operations and teleportation. 

As an application of the technique for sculpting an arbitrary motional ionic state, com- 
bined with the phenomenological approach to the analysis of the effects of noise in the 
process, we have computed the fidelity for sculpting a truncated phase state. We assumed 
realistic values for the parameters involved. Moreover, we have proposed an algorithm to 
optimize the fidelity of a sculpted state in the presence of noise. This algorithm consists 
in maximizing the fidelity of the sculpted state (subjected to noise), with regard to the 
desired state through a convenient choice of the parameters involved in the process. We 
have also shown, through the phenomenological approach, that the noise introduced by the 
Jaynes-Cummings pulse has less effect than that coming from the carrier pulse for a number 
of phonons n < 25. Thus, whether for sculpting a motional state or for any other process 
involving a small number of phonons, such as engineering a motional qubit state cq |0) + ci |1) 



for quantum computation, as performed in |^0[, the main source of errors will undoubtedly 
be the carrier pulses. 
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Appendix A 



In this appendix we show explicitly the elements 
0, 1, 2...., j, j' =t, [ and A;, k =t, |) for the carrier pulse. 



{Cn,ii {t, (f) C^_^ (t, y?)^^ = (Cn,^i {t, if) C^^^^ {t, y;)^^ = i 1 + cos{2Qt)e 

(Cn,u (t, Cl^u = (^"'iT ^^rn,u c ^ ^ ^ ~ cos(2nt)e~ 

{Cn,ii {t, Cl^^^ (t, V))^ = it, ^) CL,u = sin(2fit)e-2™'* 

/ \ / \ ^^^^ 

(C„,U (t, ^) Cl^^^ {t, ^))^ = (CnM it, ^) Cl,u it, ^))*^ = sin (2fit) e-2™'^ 
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(C„,n it, if) 4^^^ {t, ip))^ = it, ip) Cl^^^ {t, ^))*^ = ^(l- cos {2nt) e-^^^'*) , 



j) are 



at 



For the Jaynes Cummings pulse the elements ^ J Jn,jk {t, '-p) j<k' ^) 

= \ Jm,n Jn-l,ll = ^ [Cn,m " Dn,m] , 

y J = {Jm,i\ it, Jn+l,ll /j ^ [-Cn,m+1 + ^n,m+l] , 

)j = (^m,TT V') ^n,U = ^ [^n,m+l + 5n,m+l] , 

\ _ / T . . /'+ Tt ,.^\ _ r /I i 



( {t, if) (t, 



-i,TT 



where 



An,m = cos [gt {^/n - Vm)] exp -FgH {^/n - y/mf /2 
Bn,m = COS [git + \/m)] exp —Tg^t (\/n + \/m) ^ /2 
Cn,m = sin [gt {^/n - Vm)] exp -FgH {^/n - \/rnf /2 



and 



Dn,m = sin [gt {y/n + v^] exp -VgH (y^ + V^)^ /2 

Appendix B 

In this appendix we show explicitly the elements Pn,m{ti,t'2,t^) 



Pn,m{tu t,, ts) = Ai°l, [A^.y - [l - COS (2^,) e-'^""''^ 



X 
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X [An,m {t2) - Bn,m M [l + COS {2Qt,) e'^^''''''^ (1 - 5„,o) (1 - S^fi) + 
-Ail (Aai)*^-n(20^3)e-™^*3^ 

Xe-^'"' [-An,m+1 (t2) + 5n,m+l (^2)] e"^' Sill (2^1) 6-^^"'*^ (l - 5„,o) + 



xe 



[-C„,„. (t2) - A,™ (t2)] e^'^-^ sin (2f]ti) 



-2rn2(^ 



(1 - ^n,o) + 



+Ail (A(;:))%^ Sin (20^3) e---^*3^ 



xe-^'^^ [C„,™+i (t2) + Dn,m+i (t2)] [l + COS {2Qh) e-2rf^'*i] (1 - 5„_o) + 
+Ail(Ai:),)*^sin(2Qt3)e--^*3, 

[_A,„+i,,„ (t2) + (t2)] sin (2^1,) e'^™'*^ (1 - 5„,o) + 

+Ail(Aai)*^[l + cos(2Qt3)e--^^*3- ^ 

X K+i,„,+i (t2) - 5„+i,m+i (^2)] [l - COS {2Qti) e-2r^'*i 
+ASi(A(::))%^sin(2Qt3)e-^^^*3>, 

Xe'^' [Cn+l,m {t2) + Dn+l,m (^2)] [l - COS (21]ti) e'^^"'*! 

-Ail (A(:jT^[l + ^o«(2^^^3)e-^^^*3]x 

xe''^^ [C„+i,™+i (t2) + /^n+i,m+i (^2)] e-^'^i sin (20ti) e-'^"'*^ 
-Ai°)(A^,)*^[l-cos(2Qt3)e-^^^*3- ^ 

xe^'^^ [a,^ (t2) - ^n,m (^2)] e-'^i sin (2f]ti) e-^™'*^ (1 - 6^,0) + 
-Ai°)(Aai) .^sin(2Qt3)e-^^^*3>, 



+ 



xe 



-11^52 



(^2) + Dn,m+l (^2)] 1 - COS {2Vtti) 6 



-2TQ.'^ti 



+ 



-2rn2t„ 



+Ai°) {m*\[l-cos{2QU)e 

X K,m (t2) + Bn,m (^2)] [l " COS {2VLti) e-2™'*ij + 

-Ai"^) (A(:j))*^sin(2Qt3)e-^^^*3>, 

X [An,m+i {t2) + Bn,m+i (^2)] e'^^^ sin {2nh) e-2r^'*i + 



+Ai°^(ASL) .^sin(2Qt3)e-^^^^*3>, 
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;e'^^ (^2) - ^n+l,n^ (^2)] 1 + COS (2^,) 6 



-2ro2ti 



<5m,o) + 



+Ai°) (A;^^,)*^[l + cos(2^]t3)e-™^*3- 
xe-*^^ (t2) + (t2)] e^'^i sin (2f]ti) e-^™'*^ + 

-Af (A(:!))*^sin(2r]t3)e-™^'3^ 
X (t2) + (t2)] e''^' sin (21]ti) 6-2^^^'*^ + 



X 



+Ar (A(::))*i[l + cos(2fit3)e-^^^*3^ 
X (ts) + Bn+i,m+i {t2)] [l + COS (2nti) e-2™'*i 

where yji, (p2, and are the laser pulse phases for the first carrier, the Jaynes-Cummings 
and the second carrier pulses, respectively. Note that, for F = we obtain from the above 
expressions the expected results without noise. 

Appendix C 

In this appendix we show that a higher fidelity of the sculpted state, with regard to 
the desired state, does not implicate that the shape of the Wigner distribution function of 
such sculpted state must be closely to that of the desired state. We start by showing that 
for a desired state |S) there are an infinite set of sculpted state |Aa) with the same fidelity 
JF = |(Aa|S)|^, the label A standing for a continuous real parameter or a set of continuous 
real parameters to be defined below. In order to simplify our discussion we consider a desired 
vibrational state in the two-dimensional Fock space {|0) , |1)}: |E!) = (|0) + |1)) /V2. In this 
case, A reduces to a single continuous parameter since the set of states displaying the same 
fidelity JF (represented by unit vectors in the so-called Bloch sphere of M^, in analogy with 



spin-1/2 states [|T1), satisfy |Aa) = e*"^" (A |0) + e^('^i-'^o) Vl - A^ |1)), with 

J^- 1/2 



Dq = arccos ■ 



(27) 



The components =Tr aj |Aa) (Aa| of the states |Aa) in the Bloch sphere, for i = x,y,z 
and (jj referring to the Pauli pseudo-spin operators, read 
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Ty = ±2 [A^ (A^ - 1) + - 1/2)2] ^ (28) 
r. = 2A2 - 1, 

where the signal +(— ) in Vy corresponds to positive (negative) values of 0. Next, choosing 
^ = (2 + v^)/4 ?a 0.933, it follows only two states with = (r^ = 0): |Ai/2> = 
(|0) + /2 and A^g/a) = (^|0) + |1)) /2, and an infinite set of states with 7^ 0. 

All these states, with Bloch vectors lying in the cone displayed in Fig. 5, present the 
same fidehty 0.933, but different Wigner distribution functions. Figs. 6(a,b,c) display 
the Wigner functions for the states IA1/2) and ^^/2ji respectively, showing that 
although IA1/2) and ^^/2^ present the same fidelity with respect to |5), Figs. 6(b,c) 
exhibit completely different shapes compared to Fig. 6(a). It is worth mention that the 
Wigner function varies continuously when |Aa) evolves continuously from IA1/2) to K^i^ 
(going, in positive z direction, through the bottom to the top of the cone in Fig. 5, in 
clockwise or anticlockwise direction depending on the signal of Vy). 
The same analysis holds for the statistical mixtures 

px{x) = A^ |0) (0| + (1 - A^) |1) (1| + xAVT^ (e-^"^ |0) (1| + e"t> |1) (0|) , (29) 

where the parameter x e [0, 1] accounts for the purity of the density operator px{x), in a 
way that x 7^ 1 imposes thatTr {[px{K)f] < 1. For a given value of x, it follows an infinite 
set of statistical mixtures presenting the same fidelity T = (S|p;^(x)|5) when 

(b — arccos , 

From a particular choice of x = 0.9, T — 0.8, and A = 0.7, we obtain a statistical mixture 
displaying a Wigner distribution function, depicted in Fig. 7, whose shape seems to be 
closely to that in Fig. 6(a) than those in Figs. 6(b,c), although its fidelity is smaller than 
0.933. 

When considering the truncated Fock space with more than two dimension, say 
{|0) , |1) , |2)}, it follows a set of two real parameters for describing the set of states in a 
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sphere of R^, and so on. Finally, we mention that our optimization protocol described in 
Sec.V is based on the maximization of the above-defined fidelity, the overlap between the 
sculpted and the desired state, which represents a particular measure of distance between 
vectors in Hilbert space, and can not be inferred from the Wigner distribution function. 
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Figure Captions 

FIG. 1. Electronic energy level diagram of a trapped ion interacting with the laser beans 
of frequency ui and U2, where 6 = ui — uj2 — ^^o (5 ^ A), |r) (adiabatically eliminated) 
is an auxiliary electronic level which indirectly couples the levels |t) and ||), and \d) is an 
electronic level used to measure the fluorescence emission . 

FIG. 2. A quantum algorithm notation |^ for the process of sculpting arbitrary motional 



states. The complex parameters jSk and indicate rotations of the electronic states by the 
first and second C pulses, respectively, in the kth cycle and are adjusted by an appropriate 
choice of the duration and phase of the laser field following Eqs. (|5a|j5bD . U (r^, ipk) indicates 
the evolution operator for the kth JC pulse. At the end of each cycle the measurement of 
absence of fluorescence (projection on to |t)) is required for the successful accomplishment 
of the engineering process. The time proceeds from left to right, as usual. 

FIG. 3. Wigner distribution functions for (a) the initial coherent state associated with 
Ua = 0.25 and for (b) the sculpted state after the first cycle. We have used the values 
El = 25.6159 - / X 0.0379 and A = -0.3994 - / x 0.6408 x 10~^ associated with the values 
in Table 1: = 0.99, V = 0.38 and 7^ = 0.60. 

FIG. 4. Wigner distribution functions for (a) the desired phase state given in Eq. (|18|), 



(b) the mixed state sculpted in the presence of noise using the same parameters jSi and ei of 
Fig. 3(b) (JF = 0.85 and V = 0.40), and (c) the mixed state sculpted through our numerical 
optimization procedure using the set of parameters flti = 0.56, (pi = 5.48, gt2 = 0.75, 
ip2 = 1.40, nts = 1.88, and (p^ = 1.43 associated with = 0.91, V = 0.86, and 7^ = 0.64. 

FIG. 5. The cone around the axis displayed in this figure, represents an infinite set 
of unit vector states |Aa) having the same fidelity JF =(2 + v^)/4 with respect to the state 
The circles around the figure represent the contours of the unit Bloch sphere. 

FIG. 6. Wigner distribution functions for (a) the desired state |S) and the two states 
lying in the cone of Fig. 5, with = 0: (b) |Ai/2) and (c) 
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FIG. 7. Wigner distribution function for the statistical mixture px{>i) with a particular 
choice of X = 0.90, = 0.85, and A = 0.70. 
Table 

TABLE l.The probability V, fidelity JF, and rate 7^ = J^^.V^ (with ^ = 4 and C = 1/2), 
for each value of averaged excitation number Ua- The JC pulse interaction time qti and 
phase ifii in this table correspond to values which maximize the rate TZ for each ria 
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